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digital adaptive controllers 
FOR VTOL VEHICLES- -VOLUME I 
Gary L. Hartmann, Gunter Stein, and Stephen G. Pratt 

Honeywell Inc, 

Systems and Research Center 
Minneapolis, Minnesota 

SECTION 1 
INTRODUCTION 


Background 

The VTOL Approach and Landing Technology (VALT) program of the 
National Aeronautics and Space Administration is developing a technology 
base for future VTOL systems. This research includes the development 
of a navigation, guidance, and control system to permit automatic flight 
along complex four-dimensional (space and time) mission profiles from 
takeoff to landing under all weather conditions. 

One area of the overall program is developing digital control design 
procedures for VTOL aircraft and demonstrating these methods by applica- 
tion to a VALT research aircraft. This aircraft is a modified CH-47B, a 
tandem-rotor, medium-transport helicopter. Its size is representative of 
future passenger- carrying VTOL aircraft. This aircraft was previously 
used in the United States Army's Tactical Aircraft Guidance System (TAGS) 
program. * 

A large floating-point computer is part of the present on-board digital 
control system. This computer's capability permits in-flight evaluation of 


Anon. , "Tactical Aircraft Guidance System Advanced Development Pro- 
gram Flight Test Phase Report," Vols. land II, USAAMRDL TR-73-89A, B, 
Ft. Eustis, Virginia, April 1974. 


a variety of guidance and control algorithms by simply programming the 
software in FORTRAN. Control commands are transmitted from the com- 
puter to the helicopter's mechanical linkages by electrohydraulic actuators. 

As part of this research activity, Honeywell initiated work in 1977 under 
contract NASl-14921 to provide a system of digital adaptive control laws 
for flight evaluation. The control laws are self-adaptive in that they modify 
the control law on the basis of existing feedback signals. The use of air 
data sensors was specifically excluded from the control system. The control 
system was specified to be an attitude command system and was to interface 
with existing guidance algorithms. 

Candidate Adaptive Concepts 

Two different adaptive concepts were selected at the outset for evalua- 
tion using a design based on the CH-47's longitudinal axis. For maximum 
research value, one implicit method and one explicit method were selected. 

Explicit parameter estimation. -The Maximum -Likelihood Estimation 

(MLE) algorithm recently developed for the F-8 Digital Fly-by-Wire pro- 

2 

gram was selected as the explicit concept. 

This algorithm is superior to other explicit identification concepts in 
several respects. It is not subject to estimation bias in the presence of gust 
inputs which is a problem for designs based on weighted least- square s- 
output error-matching techniques. The MLE procedure differs from the 
output-error least-square technique by minimizing a residual sequence 
rather than the measurement error. By using a Kalman filter, the random 
inputs (gusts) are properly accounted for. The algorithm can also be tailored 
for real-time parameter estimation. 

^Hartmann, G. L. , et al. , "F-8C Adaptive Flight Control Laws, NASA CR- 
2880, September 1977. 
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The MLE algorithm has been implemented in a parallel- channel con- 
figuration to eliminate the need for multiple on-line iterations and eliminates 
the need for convergence tests. 

In an explicit identification approach, an algorithm is required for ad- 
justing the control law on the basis of parameter estimates. For application 
to the VALT helicopter, the control laws were structured to provide an 
attitude command system. Estimates of helicopter parameters (obtained 
on-line from MLE) perform the gain adjustment. 

Implicit model reference. -An implicit model-following design was the 

second approach investigated. Model-following systems, per se, have been 

popular for flight control, since it is often convenient to express the desired 

3 

closed-loop response of the vehicle in the form of an idealized model, 

VTOLs are no exception. As parameters in the plant vary, adjustments to 
the feedforward and feedback gains are made to maintain a response closely 
matching the model. The model used was formulated from the attitude com- 
mand response specifications. 

Document Organization 

This report is divided into 11 sections. Section 2 contains the list of 
symbols used throughout. Ground rules and requirements for the adaptive 
designs are set forth in Section 3. The helicopter model used for design is 
given in Section 4. The control structure common to both adaptive concepts 
is described in Section 5. Sections 6 and 7 contain design details for the 
MLE and the Model reference concepts, respectively. Section 8 compares 
the concepts and presents the design recommended for flight test. Section 


3 

Landau, I. D. , "A Survey of Model Reference Adaptive Techniques- -Theory 
and Application, " Automat ica, Vol. 10, p. 353, 1974. 
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9 covers the guidance interface. Section 10 summarizes results of the off' 
line processing of flight data. Conclusions and recommendations are pre- 
sented in Section 11. 
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SECTION 2 
LIST OF SYMBOLS 


Operators 

Are min f(x) 
x 

E(- ) 
exp(- ) 
ln(- ) 

In det(‘ ) 
p(x/y) 


(•) = 
A(-) 
v(. ) 


dt 


v 2 (. ) 


v p (.) 


Ml 


M 


- minimizing argument of function f(- ) 

- mathematical expectation 

- exponential 

- natural logarithm 

- natural log of determinant 

- conditional probability distribution of x given y 

- Laplace operator 

- delay operator 

- time derivative 

- increment 

- gradient vector with respect to parameter 
vector £ 

- second partial derivative matrix with respect 
to parameters £ 

- pth component of V(. ) 

- summation 

- absolute value 

T 

- quadratic form x Mx 
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Superscripts 


A 

( ) 

o 

( ) (i) 

Upper Case Symbols 

A 

b 


VZ 


D 

E q' E 0' E 
H 

I 

J 

K 

K l’ K 2' 

P 


S 

u 

V 

Y 


N 


N 


Lower Case Symbols 


a , a , a 
x’ y’ z 

c, m 


- estimated value 

- one- step predicted value 

- value for parallel channel i 

- system dynamics matrix 

- (1) discrete system input matrix 
(2) residual covariance matrix 

- measurement matrix, y due to u 

- model-following errors 

- measurement matrix 

- identity matrix 

- partial likelihood function L - l/2Z)LndetB 

- Kalman filter gains 

- control gains 

- Kalman filter covariance matrix 

- a priori parameter covariance matrix 

- Kalman filter design statistic 

- sequence of N control inputs 

- air speed 

- sequence of N measurements 


- acceleration along body axes 

- parameter vector with components c. 
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Lower Case Symbols 
(cont. ) 


g 

h 


i’- 


ii 

P 

q 

r 


t 

u 

v 

v g 

X 

y 

Greek Symbols 
Upper Case 

r 


Lower Case 
a 


0 

0 


g 


- a priori estimate of c 

- gravity 

- altitude 

- index of the minimum- L channel 

- system order 

- roll rate 

- dynamic pressure 

- (1) yaw rate 

(2) number of measurements 

- time 

- control input vectors 

- lateral velocity perturbation 

- lateral gust component 

- state vector 

- measurement vector 


- discrete system noise input matrix 


- angle of attack 

- gust angle of attack 

- angle of sideslip 

- gust angle of sideslip 



Greek Symbols 
Lower Case (cont. ) 



Ae k 


C 


r) 

e 

n 

v 

5 

P 

o, 

T 

0 

* 

CO 


- aerodynamic surface positions i = B, C, S, R 

- generic likelihood filter states 

- ij element of T 

- dummy argument for values of parameter vector 
c 

- white noise process 

- pitch attitude 

- exp (-At / t) 

- Kalman filter residuals 

- white noise process 

- air density 

- standard deviation of variable x 

- time constant 

- roll attitude 

- yaw attitude 

- natural frequency 
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SECTION 3 

DESIGN GROUND RULES AND REQUIREMENTS 


The primary application of the adaptive controller is in automatic tra- 
jectory-following. A candidate VALT flight path is shown in Figure 1. It 
is representative of the terminal portion of a VTOL's flight. It includes a 
number of heading changes and concludes with a constant-velocity descent 
followed by a decelerating descent to a hover condition. The adaptive con- 
troller will accept commands from an existing guidance law. An overall 
block diagram is shown in Figure 2, An attitude control mode was specified. 
The performance requirements for attitude control are given in Table 1. 

The adaptive control algorithm shall fit in a reasonable allocation of 
the on-board ROLM-1664 digital computer. As a design goal the adaptive 
software will operate at 10 samples per second. Sensor outputs sampled 
at 40 sps are available should a faster update rate be required. 

The sensors were limited to the present VALT sensor complement, ex- 
cluding air data. The sensors are listed in Table 2. 

The adaptive algorithm used estimates of the three Euler angles (0, <#, 
if/), the body angular rates (p, q, and r), and vertical velocity. These signals 
were supplied by existing onboard estimates that combined the various sensor 
outputs. 

Position information about the aircraft is supplied to the guidance algo- 
rithm from a barometric altimeter and a radar /laser tracking system located 
on the ground. 


The adaptive controller will operate with the unaugmented vehicle. 
Finally, the adaptive controller must not interfere with the normal mission 
of the aircraft by requiring large test signals. 


9 



East-west distance, 1000 ft 


V = 110 ft/s from® to© 

H = 750 ft from® to @ 

H = 450 ft at ©, hover at H = 50 ft at ( 5 ) 



a) Horizontal profile 



Horizontal distance, 1000 ft 


b) A ltitude profile 

Figure l.-VALT trajectory (typical). 
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Figure 2. -System functional block diagram. 


TABLE 1. -ATTITUDE COMMAND SPECIFICATIONS 


Parameter 

Requirement 

a) 

Angle commands (Q , <b. 




Rise time criterion 


Amplitude >90% of final value 
within 1.5 s. 


Overshoot criterion 


Overshoot <15% of final value. 


Settling time criterion 

Amplitude within 5% of final 
value in 5 s or less. 

b) 

V z . command: 




Rise time criterion 


Amplitude >90% of final value 
within 2 s. 


Overshoot criterion: 


Overshoot <5% of final value 
for 0 < V <10 kt. 

X 




Overshoot <0.5 V x % of final 
value for 10 < V x < 40 kt. 




Overshoot <0,5 V x % of final 
value for 10 < V x < 40 kt. 




Overshoot <20% of final value 
for V > 40 kt. 

X 

c) 

Minimal cross-axis response 
















Sensor 


Variable 

measured 


Dynamic model 


Accelerometer 

a , a , a 
x' y' z 

1 

1 

1. 005s + 1 

Angular rate gyro (3) 

P, q, r 




s + 2(0. 65) (100)s + (100)" 

Vertical gyro: 


1 


Pitch 

Roll 


0 

0 


Gyromagnetic compass 


* 


door 


s 2 + 2(0. 65) (100)s + (100) 2 


Instantaneous vertical 
speed indicator (IVSI) 


TS + 1 


Barometric altimeter 


0. 07s '+ 1 


Calibrated airspeed 


0.025s + 1 


a 


Location with respect to c.g. 


x(in. ) 
(m) 

y(in. ) 
(m) 

z(in. ) 
(m) 

Dynamic range 

10 

(0.254) 

0 

(0) 

-20 

(-0.508) 

±0. 5g (±4. 9 m/s 2 ) 

-110 

(-2.8) 

20 

(0.508) 

30 

(0.762) 

±30° (±0.52 rad) 

-110 

(2.8) 

20 

(0.508) 

30 

(0.762) 

±30° (±0.52 rad) 
±45° (±0.78 rad) 

Not specified 

0 to 360° (0 to 6.28 rad) 

Not specified 

±200 kt (±103 m/s) 

/ 

Pressure ports at 

0 to 5000 ft 
(0 to 1524 m) 

280 
(7. 11) 

40 

(1.016) 

-10 

(-0.254) 


Aircraft nose 


50 to 200 kt 

(25.74 m/s to 103 m/s) 

















































SECTION 4 

MODELS AND SIMULATION 


CH-47 Model 


The helicopter model used for analysis and design consists of linear 
differential equations and a corresponding set of stability derivatives. 


The tables of stability derivatives are taken from Appendix E (pp. 85- 
4 

93) of ref. 4. They were used to define 99 flight conditions. These are 
numbered consecutively such that each flight condition corresponds to a 
column of Table E2. The state vector, x, and the control vector, 5, for 
the model are defined below: 


State Vector 
u 

v 

w 

P 

q 

r 

0 

e 

* 


Control Vector 

6-d differential collective 
Jd 

6^ collective 
6 S cyclic 

6^ differential cyclic 


^Ostroff, A. J. , Downing, D. R. , and Root, W. J. , "A Technique Using a 
Nonlinear Helicopter Model for Determining Trims and Derivatives, ” 
NASA TN D-815 9. 
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For control design purposes, the model was decoupled into separate 
longitudinal and a lateral-directional models. For the pitch-axis model, 
the u and v velocities were rotated to earth-referenced coordinates, x and z. 
This was done to permit direct evaluation of guidance command responses. 
The variables, control inputs, and guidance commands are given in Table 3. 


TABLE 3. -CH-47 VARIABLES 


Parameter 

Pitch model 

Lateral-directional 

model 


V x 

V 

Variable 

V z 

P 


q 

r 


e 

0 

Control 

6 b 

6 S 


6 c 

6 r 

Command 

e 

0 


V z 

* 


Each of the uncoupled models is represented by a linear model in the 
usual form 

x = Ax + B6 

Thirty-three trim points (flight conditions) were selected for design and 
analysis purposes. These 33 flight conditions cover the range of the CH-47' s 
operating envelope shown in Figure 3. 
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Figure 3. -CH-47 flight conditions. 


Actuator and Rotor Models 

The rotor and actuator model used for design is shown in Figure 4. This 
configuration evolved from the CH-47 TAGS systems as modified for the 
NASA Langley VALT program. The gain was introduced to convert any 
stick deflection into standard CH-47 stick inches compatible with the set of 
stability derivatives. 


Simulations 

A digital simulation of the CH-47 was developed for use at Honeywell 
in designing the adaptive controllers. It operates with a specified sample 
rate and integrates the differential equations to compute CH-47 time histories 
in response to commands and/or disturbances. The simulation uses constant 
parameters representing a fixed flight condition. The simulation has the 
option of representing either longitudinal, lateral-directional, or coupled 
longitudinal- lateral dynamics. 

A random-number generator is used to represent sensor noise, and 
random sequences shaped by 1-second first-order lag filters are used to 
represent forward and vertical gusts. 
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The simulation includes the linear model for the rotor dynamics from 
Figure 4. Rate limiting of the actuator was included, since this nonlinearity 
tends to affect the stability of command responses and thus is important. 

The simulation used for design verification was the nonlinear six-degree- 
of-freedom VALT simulation at NASA Langley Research Center. It operates 
in real time on a Control Data CYBER computer and has the capability of 
interfacing with the actual flight computer. 

The control laws designed for this program are described in the next 
sect ion. 
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SECTION 5 

CONTROL LAW STRUCTURE 


A design procedure based on state- variable feedback was used to design 
the attitude command system. The design procedure uses an algebraic 
model- following method and includes integral control to eliminate steady - 
state errors and improve decoupling. Uncoupled controllers were designed 
for the pitch and lateral-directional axes. Each design is discussed below. 

Longitudinal Control 

The specified 6 attitude response is essentially a second-order response. 
Consideration of the desired speed of response leads to a requirement for 
to > 1. 3 rad/s and f > 0.5 in a second-order response. The closed-loop 
pole corresponding to the integral was placed at -0. 8, which is sufficiently 
fast so that the response is dominated by the complex pair of poles. If the 
integrator pole is moved further left, rate- limiting problems are encoun- 
tered. A feedforward gain was used to add a zero at -0. 8, cancelling the 
closed- loop integrator pole. Thus, the desired transfer function for a 6 
command is: 

Q_ _ 6 (s + 0. 8) 

6 c (s 2 + 3s + 6) (s + 0. 8) 

Similarly, the first-order V response has a desired transfer function: 

z 

V 2 (s + 1) 

Z _ ■ 

V " (s + 2) (s + 1) 
z 

c 
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where the feedforward is again used to cancel the closed- loop integrator 
pole. 


The addition of integral control adds two additional states. Thus, the 
pitch-axis design involves six states and two command inputs: 

X T = (V x , V z , q, Q, I v I 2 ) 

where 



and 

uT * V 


The desired transfer functions can be put in state form: 


x = A x + 13 u 
m m 


where x and u are as defined above. This yields: 


- 

- 

- 

- 

- 

- 


0 

0 

0 

-3 

0 

0 

2 

0 


2 

0 

0 

0 

-4 

-9 

0 

6 

B 

0 

6 

0 

0 

1 

0 

0 

0 

m 

0 

0 

0 

0 

0 

-1 

0 

0 


1 

0 

0 

- 0 . 

8 0 

0 

0 

0 


0 

0.8 
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to realize the desired 9 and V transfer functions. (The first row of A 

z m 

has been ignored, since V is not being controlled. ) 

X 

The vehicle model has the form: 
x = Ax + B5 

where the actuator commands (6) are 

6 = K x + K U 
u x u 

The control design task is to specify feedback gains, K , and feedforward 
gains, K u . 

The block diagram is shown in Figure 5. 



Figure 5. -Pitch-axis block diagram. 
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The feedback gains were computed algebraically from 


K x -H ' 1 [ A m- A ] < ! > 

where B is a 2 x 2 submatrix containing rows two and three of B, and rows 
two and three of A m - A are used (indicated by overbar). 

Similarly, 

K u"[ g ]‘ lB m ■ < 2 > 

The 12 elements of the gain matrix were computed at 33 flight conditions 
using Equations (1) and (2). The values are plotted in Appendix A. An 
analysis of the gains indicated they could be made simple functions of V . 
These schedules are given in Table 4 and also plotted in Appendix A. For 
convenience, the scheduling parameter was defined as forward velocity 
normalized by its maximum value (v). Thus, -0.25 £ v £ 1.0. 

Transients responses for 0 and V commands based on these gain 

z 

schedules satisfy the specification and exhibit very little cross coupling. 
Figures 6 and 7 show a representative response and the specification enve- 
lope for a V and a 0 step command, respectively, 
z 

Lateral-Directional Axes 

A lateral-directional controller was designed for the CH-47 to accept 
0 and guidance commands. The feedback and feedforward gains were 
computed to provide the desired command responses and to minimize the 
cross-axis response. The controller incorporates integral control for each 
command and uses measurements of p, r, 0, and i// in the feedback. The 
controller outputs are 6 (cyclic) and 6„ (differential cyclic) commands. 

S IX 
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The lateral-directional controller also has the functional block diagram 
shown in Figure 8, where 

T 

6 = (5 , 6 D ) are the actuator commands 

S ri 

T 

u = (0 , i// ) are command inputs from guidance 
c c 

T 

X = (v, p, r, 0, 1 p, I 3 , I 4 ) 

and two command error integrators have been appended to the vehicle 
states: 



h = ^c'^ 


TABLE 4. -PITCH-AXIS GAIN SCHEDULE 


K x 

element 

Function 

(1.1) 

0 

(1. 2) 

-0. 25 v v < 0. 5 


-0. 125 v * 0.5 

(1.3) 

-7.6 + 1.5 v 

(1.4) 

-27+6 v 

(1.5) 

0. 14 v v < 0. 5 


0. 07 vS0.5 

(1, 6) 

17.5 -5 v 

(2, 1) 

0 

(2,2) 

0. 3 -0. 12 v 

(2,3) 

-l.Ov v<0.5 


-0. 5 v s 0. 5 

(2,4) 

-18 v 

(2,5) 

-0. 23 + 0. 08 v 

(2,6) 

2 v v < 0. 4 


1. 04 - 0. 6 v vi 0.4 
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Normalized V . Normalized 










The desired gain values were determined in a manner analogous to the 
design of the pitch axis. The desired transfer function for 0 and t// com- 
mands is identical to the 6 transfer function. A model with system matrices 
A m and B m was defined that provides uncoupled <t> and responses that 
satisfy the transient response specification. The gains were determined by 
the relations: 

K x * N'lv 1 ! 

K u = [s^KJ 


As the CH-47 parameters (A, B) change with flight condition, the de- 
sired gain will vary. Appendix B presents a pLot of the desired gains for 
33 flight conditions covering the CH-47' s flight enveLope. The gains are 
defined as: 


V 


K 11 

K 12 

K 13 

K 14 

K 15 

K 16 

k i 7 ' 

1 

o> 

w 

1 


K 21 

K 22 

K 23 

K 24 

K 25 

K 26 

K 27_ 


v 

P 

r 

0 

Io 


where Ig 


and I . 
4 


are the command error integrators. 


It is worth noting that there is no strong dependence on forward velocity 
as there was in the pitch axis. The variation of gains over the 33 flight con- 
ditions is also rather limited. The major trend in the data seems to be a 
variation of gains with vertical velocity for forward speeds greater than 
41. 1 m/s (135 ft/s). This is probably due to similar variations in elements 
of the B matrix. 
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As a first design, the 14 gains were held constant over the flight 
envelope. (Gains on sideslip were set to zero.) The responses to 0 and 
\f/ commands were evaluated over the flight envelope using a coupled pitch 
and lateral simulation. The 0 and i p responses meet specifications and 
exhibit minimal cross-axis response even at forward airspeeds greater 
than 41 m/s. In addition, no interaction was observed between the pitch 
and lateral-directional controller. Representative 0 and time histories 
are contained in Appendix C. A representative response is shown in Fig- 
ure 9. 

It is concluded that fixed lateral gains are completely satisfactory 
and there is no advantage to adjusting them. Thus, the system recommended 
for flight test consists of an adaptive pitch controller and a fixed-gain 
lateral-directional controller. 

The next two sections describe details of each of the adaptive algorithms 
evaluated. 
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SECTION 6 

THE MAXIMUM -LIKELIHOOD A LGORITHM 


This section presents the design of an explicit parameter estimation 
algorithm. It was combined with algebraic gain computation to serve as 
one of the candidate adaptive concepts. 

First, a review of the MLE Algorithm is presented, then the design 
issues are discussed. This is followed by a summary of the simulator per- 
formance of the algorithm. 

Overview of Algorithm 

Parameter uncertainties in aircraft generally take the form of unknown 
parameters in an otherwise known model structure (i. e. , coefficients of 
linearized equations of motions). Their range of uncertainty is largely due 
to widely varying flight environments (dynamic pressure, velocity, angle of 
attack) and configuration variations (center of gravity, fuel, payload). How- 
ever, the individual coefficients are strongly interrelated and only a few 
must be known accurately for control law design. The explicit identification 
problem then is: 

Given a continuous plant 

x = A (c)x + B(c)u + T(c)5 

with discrete measurements 
Y. = H(c)x. + D(c)u. + n ^ 
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find the unknown parameter vector, c. The 5 is a white noise process 
driving gust states and the p. are sensor measurement noise. 

There is a variety of explicit identification methods. If the inputs and 
outputs of a plant are known, a least- squares equation error is appropriate. 

If the outputs are noisy, then an output error method will provide unbiased 
estimates. Finally if the plant has stochastic inputs and noisy measurements, 
a Kalman filter approach provides unbiased estimates. 

If the parameters are modeled as states and the resulting filter is non- 
linear, approximations can be implemented (see ref. 5 for an aircraft 

5 

example). If one is interested in the first and second moments, then a 
maximum likelihood procedure can be used. This approach was attractive 
because it provides unbiased estimates (in the limit) with random inputs and 
noisy measurements. 

The general process of maximum -likelihood identification is to find 
parameter estimates which maximize the a posteriori probability distribu- 
tion for the observed outputs conditioned on the unknowns and the measured 
inputs; i. e. : 


c = Arg [max p(y r y 2 , . . . y n | c = Q; u Q , ' u r . . . u nl )} 

When the unknowns are constant and the plant dynamics are linear, this 

6 

maximization problem leads to the solution shown in Figure 10. The 


5 

Gelb, et al. , Applied Optical Estimation, The MIT Press. Cambridge. 
Mass., 1977. 

O 

Balakrishnan, A.V. , "identification and Adaptive Control: An Application 
to Flight Control Systems," J. Optimization Theory and Applications, Vol. 
9, No. 3, March 1972. 
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solution consists of a Kalman filter designed for the true system structure 
but with parameters equal to an estimate, c = £. The filter tracks the true 

A 

system outputs and generates a residual sequence {v^ = k = 1» 2, 

This sequence is accumulated into a likelihood function: 

L (C, N) ^ - In (y 1 , y 2 , • • • Y n I 9 = Cs u 0 » • • • u n _ 1 ) 

which is then minimized over the parameter estimate, £. At first glance, 
this solution appears ideal for onboard applications. The Kalman filter is 
relatively simple and runs recursively, processing data samples as they 
appear. The same is true for the likelihood accumulation operation. The 
difficulty, of course, is the last step of the solution: the minimization. This 
requires repeated or parallel processing through the data and adds signifi- 
cantly to computational complexity. Two algorithms were considered for 
the likelihood minimization operation: 1) iterative Newton- Raphson cal- 
culations, and 2) parallel noniterative calculations. 


Inputs 



Figure 10. -Maximum- likelihood estimation. 
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The first algorithm begins by collecting a sequence of input/output data. 
The Kalman filter is then run with parameters c equal to a priori estimates, 
c = £° = c q . This generates the likelihood function, L. At the same time, 
a set of sensitivity equations is processed which permits calculation of first- 
and second-partial derivatives, vLand v L (some form of approximation 
is usually used to simplify equations for the latter). These derivatives are 
then used to obtain a new parameter estimate using a standard Newton- 
Raph son formula. The filtering, accumulation, and Newton- Raphson opera- 
tions are performed repeatedly for the same data set until convergence is 
achieved. The data set itself is usually kept current by a "sliding window" 

process. This algorithm has been implemented successfully for various 

6 

postflight data processing applications. 

The second algorithm replaces iterative calculations with parallel ones. 
The sequence of input/output observations is sent simultaneously to M Kal- 
man filter channels, each with its own sensitivity calculations and likelihood 
accumulations. The channels are distinguished by their assumed parameter 
values. Each one operates with a different parameter estimate, c = 
and hence computes, L^, VL^, and at a different fixed point in param- 

eter space. The likelihood functions at these points are then compared to 
find the approximate minimum point, and a single Newton- Raphson step 
is taken from there to estimate the true extremum. As in Algorithm 1, 
some procedure must generally be added to keep the processed data current. 
This is done by highpass operations which provide exponential de-weighting of 

past accumulated data samples. General parallel estimation structures have 

7 8 

been suggested in several references. ' 


7 Stein, G. and Saridis, G. N. , "A Parameter Adaptive Control Technique, " 
Automatica, Vol. 5, pp. 731-739, 1969. 

8 

Lainiotis, D. G. , " Optimal Adaptive Estimation: Structure and Parameter 
Adaptation," IEEE Trans. Auto. Control, Vol. AC-14, No. 2, April 1959. 
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As criteria for onboard application, the recursiveness and fixed struc- 
ture of Algorithm 2 make it ideally suited for real-time operation. 

Filter, Sensitivity, and Likelihood Equations 

The basic data processing equations which must be resolved by each 
filter channel are summarized in this subsection. They are stated in terms 
of general symbols corresponding to the following discretized plant equa- 
tions: 


x k+i + Ax k + Bu k + r,i k 
y k = Hx k + Du k + N5 k 

The matrices A, B, T, H, D, and N should all be thought of as dependent 
functions of the parameter vector c. Then the channel equations are: 


Filter equations: 

x k+i = A5 k + Bu k 

v k+i ■ y k+i - Hx k+i - Dv k+i 

x k+l " x k+l + K k+1 v k+l 

Sensitivity equations for each component c^ of c (derived by dif- 
ferentiating with respect to c ): 


Vk+i * A Vk + ( V )x k + (7 P B>u k 


Vk+1 ■ - H7 p 5 k + 1 - (7 p H)x k + i - (7 p D>u k + i 

7 p x k+l ~ 7 p x k+l + ^k+1 7 p v k+l + * 7 p^k+l * J k+l 
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2 

Filter gains: 

P k + 1 

= a p k a + rr T 

B k+1 

= (H P k+1 H T + N N T ) 

K k + 1 

■ f k + l h G W 1 

P k + 1 

= P k+1 ■ K k+1 H P k+1 

Like lihood 

accumulation: 

L k + i 

■ ivii i 2 w 1 

7L k+i 

= MVL k+ W k+1 T B k+1 _1 v k+1 


+ 1/2 Trace (\ +1 ' , k+1 T - B k+1 ) 7 ®k+l ' 

^ Vl = ^ 2L k + W k + 1 T Bfell' 1 77 k + l 


with p = exp (-At/r) for exponential deweighting of past data. The choice of 
r is discussed in this section. 


These equations warrant two explanatory comments. First, there are 
no sensitivity equations for the filter gains. This is because the matrix VK 
was computed by numerical differentiation throughout the design program; 
i. e. : 


V K = [K(c + \ e ) - K(c)]X 
P - P P ~ P 


where e is a unit vector in the ith coordinate direction and ^ was chosen 
P P 

small compared to the range of c . 

P 
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2 

Second, since v L includes approximations, all second-partial deriva- 
tives and products of derivatives have been ignored. This is a common 
approximation for so-called modified Newton-Raphson procedures and has 
the important advantage of eliminating second -partial derivative sensitivity 
equations. It is also common practice to eliminate the last term of the VL 
equation and to replace the filter gain equations with their steady- state 
solutions. 


Parameterization 

The CH-47 model was parameterized for several reasons: 

• Some "structure" must be imposed on the model if the identi- 
fication is to work at all. If each of the elements of the A 

and B matrix are treated as independent variables, one simply 
cannot learn much by processing 5 to 10 seconds of data. 

• Real-time computing constraints demand that only a limited 
number of parameters be identified. With a software structure 
like PCMLE, the computer time and memory required grow 
linearly with the number of parameters we try to estimate. 

• From a control viewpoint, only a few "dominant" characteris- 
tics are important for adjusting the control law. For the 
CH-47, as shown in Section 5, the v parameter turns out to be 
the most important factor for this function. 

The CH-47 pitch axis model was parameterized by expressing all the 
A and B elements in terms of the v parameter plus individual perturbations 
to this function. The initial uncertainty was estimated from the range of the 
parameters over the flight envelope. Identif iability then indicates how much 
reduction in the initial uncertainty can be expected by processing practical 
amounts of sensor data from a physically realistic noisy environment. 
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Plots of the longitudinal model coefficients are contained in Appendix D 
for 33 flight conditions that completely cover the CH-47's operating envelope. 
The parameters are plotted against forward airspeed (-20. 5 to +82. 3 m/s) 
for three values of vertical velocity (-10. 0 m/s, 0, and +10. 3 m/s). Each 
of these 18 model parameters was characterized in a manner analogous to the 
F-8 Adaptive Study. Parameters showing a strong dependence on normalized 
forward airspeed (v) were made a function of v plus a perturbation parameter 
c.. If the parameter didn't show a v dependence it was defined as a constant 
(which could be zero) plus the perturbation quantity. This method of param- 
eterization is used to reduce the initial parameter uncertainty. The param- 
eterization being used is given in Table 5 and is also plotted in Appendix D. 

A model based on this parameterization was used in the following identifia- 
bility analyses. 


Identifiability Analysis 

The design issues are primarily resolved with an "identifiability 

analysis. " A linear system identification problem is formulated to estimate 

a parameterized model from closed- loop input/output records generated by 

9 

a simulation. As shown in the literature, the theoretical accuracies attain- 
able are bounded from below by the "Fisher information matrix. " Thus, this 
matrix is computed for various operational situations and its accuracy bounds 
are used to structure the identifier. 

The model shown in Figure 11 was used for the pitch-axis identifiability 
analysis. This model includes first-order gusts models for forward and 
vertical gusts. The bandwidth of each gust process is one rad/s. 


^Tse, E. , "information Matrix and Local Identifiability of Parameters," 
JACC, Columbus, Ohio, 1973. 
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Function 


RMS c 


a ll 

= 

-0. 018 - 0.034 v 

+ C 2 

0. 0045 

a 12 

= 

.2 

0. 048 v + Cg 


0. 009 

a 13 

= 

2. 8 + c 4 


0.4 

a 14 

= 

-32.2 + (14 + c 5 ) 

-2 

V 

3. 0 

a 21 

= 

(0.18 + c 6 )M^ 

^ s,nv 

0. 04 

a 22 

- 

-0. 5 + Cr^ 


0. 25 

a 23 

= 

C 8 


1. 6 

a 24 

= 

-180 V + Cg 


10. 0 

a 31 

= 

c 10 


0. 02 

a 32 

= 

0. 02 v + c 


0, 006 

a 33 

= 

-1. 5 + c 12 


0. 25 

a 34 

= 

4 + C 13 


0. 5 

b U 

= 

0. 12 + c 14 


0. 06 

b 12 

= 

C 15 


0.3 

b 21 

= 

0. 70 v + c^g 


0. 25 

b 22 

= 

-7. 8 - 3v + c 17 


0. 8 

b 3 1 

= 

0, 35 4-0.12 v + 

0. 05 

b 32 

S 

0.24 v + Cjg 


0. 08 









Figure 11. -CH-47 model for identif iability analysis. 


For the majority of the analysis, the sensor rms noise magnitudes 
were selected as: 

CTy = 1 . 3 m / s 
z 

(7 = 0. 0026 rad /s 

Oq - 0. 017 rad 

The CH-47 was modeled with v as a dominant parameter plus smaller 
perturbation parameters yielding 19 parameters. Adding two parameters 
for X and Z gust intensity plus the six initial conditions on the state results 
in a total of 27 longitudinal parameters. 
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The identification problem was then formulated to estimate these param- 
eters from closed-loop input (6,,, 6 ) and output V , q, 0) time histories. 

ij c z 

For the majority of runs, 5 seconds of data were used. 

Qualitative characteristics of the identifiability analysis are given in 
Table 6. The results basically show that under most conditions about five 
to nine of the perturbation parameters can be identified (parameter uncer- 
tainty reduced by at least a factor of two). The quality of the identification 
does not vary much with flight condition. The gust levels can be estimated 
only when they dominate command and sensor noise effects. (This was 
observed on F-8 also. ) 

Assuming reduced sensor noise does not dramatically change the results, 
note that for many of the model parameters, knowledge of v yields a good 
estimate. The perturbation parameters with very small uncertainties are 
not important. Experiments with estimating a subset of seven parameters 
have shown little degradation in their accuracy due' to the parameters not 
being estimated (also observed on F-8), 

In our initial identifiability analyses, most of the data used both A0 and 

AV commands and predicted accuracy results were good. However, when 
z 

AV commands only are used, the predictions get substantially poorer, par- 
z 

ticularly when reduced -parameter identification is attempted. Results for 

this AV -only case are presented below, 
z 

Figure 12 shows the reduction in the v uncertainty using 5 seconds of 

data containing AV commands only. The reduction is shown as a function 

of the parameter set identified. Even with all the parameters estimated, 

the theoretical bound is ct =0. 15, which is not too good. Our baseline 

parameter set has a = 0. 3, and a better four -parameter set for V com- 

z 

mands would give o = 0. 17, which is nearly as good as estimating all the 
parameters. 
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TABLE 6. -SUMMARY OF PITCH-AXIS IDENTIFIABILITY 


Flight 

condition 

Command 

RMS 

gusts 

Sensor 

noise 

Number of 
paramested 
estimated 

Parameters identified 

Comments 

55 

A 0 = 10° 

0 

No 

All 

V ' a 3T a 33' b 31' b 32 

Command is doublet with 
4-sec period. 

55 

AV = 10 ft/s 
z 

0 

No 

All 

v, a 2 2 » < a 32 ) » b 31* b 32 

Command is doublet with 
4-sec period. 

55 

A 0 = 10° 

AV * 10 ft/s 
z 

u g - 1 ft/s 
W g = 1 ft/s 

Yes 

All 

v, (a 22 ). a 31 , a 32 , a 33 , b 31 ,b 32 

Commands push gust estimate 
to lower than actual value. 

55 

A 0 = 10° 

AV = 10 ft/s 
z 

L_ . 

U g = 1 ft/s 

w = 1 ft/s 
g 

Yes 

7 Identified 
above 

v ' a 31* a 32* a 33' b 31* b 32 

a 22 and a 3 j degraded when 
subset estimated; a 3 j still 
improves. 

55 

Ae = io° 

AV = 10 ft/s 
z 

u = 1ft /s 
8 . 
w = 1 ft/s 
g 

Yes 

All 

v, (a 33 ) # b 31 , b 32 

Only 1 s of data used. 
Response on v ,b 31# b 32 

is very fast. 

55 

A 0 = 10° 

AV = 10 ft/s 
z 

u g = 1 ft/s 

w = 1 ft/s 
g 

No 

All 

v, (a 22 K a 31 , a 32 , a 33 >b 31 ,b 32 

Gust estimates go even closer 
to zero than row 3. 

55 

A 0 = 1° 

AV *=1 ft/s 
z 

Ug = 1 ft/s 
W g = 1 ft/s 

No 

All 

v, (a 31 ) # (a 33 ), b 31 , b 32 

A z = 0. 06g rms 

0 = 0.75° rms 
V z = 0. 8 ft/s rms 

55 

None 

U = 1 ft/s 

w = 1 ft/s 
g 

Yes 

All 

(v). a 31 

A = 0, 047 g rms 
z 6 

0 = 0. 025° rms 

V = 0. 55 ft/s rms 
z 

55 

Random A 0 
Random AV 

z 

o 

No 

All 

(v) # b 31 , (b 32 ) 

A^ = 0. 05g rms 

0 = 0, 083° rms 
V z = 0. 17 ft/s rms 

Note lower bandwidth of V and 
z 

0 loops results in less excitation 
with random inputs. 

55 

i 

Amplitude 
doubled from 
above 

0 

No 

All 

v, (a 33 ), b 31 , b 32 

More low-frequency test signal 
would be better. 

5 

A 0 = 10° 

AV * 10 ft/s 
z 

u = 1 ft/s 
w g = 1 ft/s 

Yes 

All 

v, (a 22 ) f a 31 , a 32 , (a 33 ) 
(a 34 ), (b 22 ), b 31 , b 32 

Doublet commands compare with 
row 3 (Fc 55) 

93 

A0 = 10° 

AV = 10 ft/s 
z 

U = 1 ft/s 1 
2 * 
w = 1 ft/s I 
g 

Yes 

All 

v# a 22* a 31* * a 32 J ' a 33 
< a 34> # ^ a 22^* b 31* b 32 

Doublet commands compare with 
row 3 (Fc 55). 

47 

A 0 = 10° 

No 

No 

All 

v, a 31 , (a 33 ), < a 34 ), b 31 

Doublet commands 

47 

AV = 10 ft/s 
z 

No 

No 

All 

v, a 22 , (a 23 ), a 32 , b 21 ,b 22 ,b 32 

Doublet commands 

47 

A 0 = 10° 

AV = 10 ft/s 
z 

u = 1 ft/s 
6 

w = 1 ft/s 
g 

Yes 

All 

v, (a 22 ), a 31 , a 32 , a 33 , (a 34 ) , 
<b 22 )# b 31* b 32 

Doublet commands 

' 47 

AO - 10° 

AV = 10 ft/s 
z 

u g - in ft/s 

w =10 ft/s 
g 

Yes 

All 

Same as above plus u and w 
gust Intensity 

Doublet commands 


( ) Indicates Improvement In Initial uncertainty 
but by less than a factor of two. 
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Initial 



Fc 51 (80 kt) 

M/ z =10 ft/s doublet 

5 s of data 

No gusts or sensor noise 


Parameter sets 
1 = v 

4 = v, a 33 , ^ 21 ' ^32 baseline set ^ 

4' = v, a 32 / ^31' ^32 


5 

13 


— v, a 32' a 33' ^31 # ^32 

= v, a 12 , a 22 , a 31 , q 33 , a 34 , b 21 , b 22 , b 31 , b 32 , u g , w g 


Figure 12. -Theoretical 1- sigma accuracy in v as function of 
parameter set estimated. 


Subsequent analysis of the likelihood function indicated that the pitch rate 

residuals are heavily weighted. Much of the identification is therefore based 

on the pitching moment response. For AV responses, the coupling term 

z 

a 32 is important. Figure 13 shows the influence of this parameter on v 
accuracy for the baseline set of four parameters. (This confirms the dif- 
ference between set 4 and 4' in Figure 12). 


On the basis of Figure 12, several runs were made with the PCMLE 
algorithm modified to estimated a^ rather than a^g. There was no improve- 
ment in the estimate for the V x = 135 ft/s flight condition, although (theore- 
tically) the 1-sigma accuracy bound is smaller. If the A and B matrices cor- 
respond to the parameterization, then v is correctly estimated. A number 
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of experiments with varying elements of the A and B matrices lead to the 
conclusion that errors in parameters not estimated shift the minimum of the 
likelihood function from the correct value of v. Furthermore, this effect is 
not caused by one or two parameters but is the effect of all the parameters 
contributing to the uncertainty of v. 

.5 


CT - 
V 

0 .0006 .0024 .006 .01 

°a 

a 32 

Figure 13. -Uncertainty in v as a function of uncertainty 
in a 32 model coefficient (four-parameter 
identification, V z commands). 

For comparison, Figure 14 shows the reduction in v uncertainty with A0 
doublets (no gusts or sensor noise). The baseline set of four has or = 0. 12 
compared with cr = 0.075 for estimating all the parameters. Set 4', which 
looked better with V z commands, does not perform as well as the baseline 
set for A 0 commands. 

Finally, Figure 15 shows that the theoretical accuracy level is relatively 
constant over the flight envelope. 

Analysis of Likelihood Function Contours 

An important aspect of MLE parameter estimation concerns the shape 
of the likelihood function away from the correct parameter values. This 
significantly influences the region of parameter space that can be estimated 
from any given channel. 
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o 1 4 5 7 all 


No. of parameters identified 

Parameter sets 
1 = v 

4 = v, a 33 , b 31 , b 32 
4' = v, a 32 ' ^31' ^32 

5 = v, 3j2' a 33' ^31' ^32 

7 = v, a 31 , a 32 , a 33 / a 34» b 3p ^32 

Figure 14. -Theoretical 1- sigma accuracy in v as function of 
parameter set measured. 



Figure 15. -Theoretical uncertainty level for estimating all 
parameters as a function of flight condition. 
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Representative contours produced by the initial filter designs are shown 
in Figure 16. The major problem is near hover. Here a very sharp notch 
results because the filter has the forward and vertical gust states decoupled 
from q and 6 states at hover. This problem was remedied by redesigning the 
filter. Figure 17 shows PCMLE estimates from simulation data based on our 
v parameterization of the CH-47 pitch axis. Note the minimum occurs 
at the correct parameter value, the shapes are approximately quadratic, and 
the gradients computed away from the true parameter point to the correct 
answer (parameter value corresponding to minimum). 

By plotting likelihood functions, it was determined that v could not be 
estimated accurately by minimizing the likelihood function with this single 
parameter. This effect is shown in Figure 18 which illustrates minimizing 
only along the v direction, assuming the perturbation parameters (c.) are 
zero gives an erroneous v estimate. Thus, it was found to accurately esti- 
mate v required estimating three additional parameters. This is discussed 
further in the next section. 

Identifier Design for CH-47 

An identifier was developed for the pitch axis in view of the identifiability 
results and the control requirements for gain adjustment. The identifier 
uses the parallel minimization algorithm previously discussed. It is based 
on a reduced parameter set and identifies four parameters -- v plus pertur- 
bations on B(3, 1), B(3, 2), and A (3, 3) (recall state equations of vehicle). 
Four state, constant-gain Kalman filters are used. Three filters are used 
in the baseline configuration to cover the operating envelope. 

The major design issues include: 

• The identification model 

• Channel selection 
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• Kalman filter design 

• Adaptation to noise statistics 

• Likelihood filters for deweighting old data. 


Each of these items will be briefly discussed. 


Contours of constant- 



True minimum at: 
v = 1 

c i = c * 


Figure 18. -Multiparameter likelihood contour. 


Identification model. -Pitch-axis dynamics are modeled as: 
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The parameters Y^ and Y 22 represent gust disturbances as defined 
below. 

The Dryden gust spectrum 

■ "X v g + a XT T ' 
was approximated by 

V =- cr -7- n 
g L 

A gust model was used in the x and z axes. 

The bandwidth of the gust filter was set at 1 second; thus: 

Y n 0 Xg 


Y 22 = \/F a z 


g 


The Ygg term was added to improve the shape of the likelihood functions as 
described previously: 

Y 33 c a~ (0 - 02 a x +0 ' 01 ' J z > 

33 g g 

Measurements were modeled as 
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The sensor noises are modeled as independent and identically distributed 

random variables. Their magnitudes, ct v , a a , and gq, as well as the dis- 

v z 4 

turbance magnitudes, o x ^, a z * were treated as design parameters. 

The trim values of the states and controls were eliminated from the 
identification problem by highpassing all the data with a second-order filter. 
The break frequency was set at 2 rad/s. 

Channel selection. -The problem of channel selection is to choose both 
the number and location of points in parameter space. Some insight into this 
problem is provided by the shape of the likelihood function away from the 
true parameter value. The shape and curvature indicate how well interpola- 
tion can be done from each channel. Some experimentation with the simulation 
is also useful in selecting the channels. 

The baseline design used three channels spaced at v = 0, 0. 5, and 1. 0, 
as shown in Table 7. 


TABLE 7. -NOMINAL CHANNEL LOCATION 




Parameters 

in Kalman Filter 

Channe 1 

c-^v) 

°2 < ^ a 33 > 

C 3^ Ab 31* 

c 4^ b 33^ 






1 

1.0 

-0. 045 

-0. 05 

r— 

o 

o 

1 

2 

. 5 

0 

0 

0 

3 

0 

0 

0 

0 


Kalman filter design. -According to Algorithm 2, a fixed set of Kalman 
filter and sensitivity equations must be operated at each set of parameter 
values. The design of these filters is straightforward, since steady-state 
gains can be used. Hence, the gains are solved off-line for their steady-state 
values. Stored values of K and VK are used in the on-line identifier. 
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Adaptation to proportional noise statistics. -We noted in the filter design 
section that fixed statistics were used to compute filter gains. This is de- 
sirable because it generates an invariant set of gains for each channel which 
can be computed off-line and stored for on-line use. Invariant gains are 
actually obtained under slightly less restrictive circumstances, namely, 
when the disturbance and noise statistics remain in constant relationship to 
one another. This means that a filter for disturbances statistics and sensor 
noise statistics remains unchanged when both statistics are scaled by the 
same factor. Asa result, the total identification algorithm designed for one 
set of statistics can still be valid when those statistics are scaled up or 
down, provided only that the likelihood functions are also scaled. This was 
done adaptively as described by Hartmann, et al. Assume that all statistics 
are known to within a constant scale factor, cr. This parameter is estimated 
and used to scale the likelihood functions before a comparison is made to 
select the minimum channel. 

Likelihood filters. -Highpass filtering of the accumulated likelihoods in 
each channel keeps the accumulations current by exponentially deweighting 
past data. The rate at which deweighting occurs (or the choice of time con- 
stants for the highpass filter) is determined by two conflicting requirements; 

1) A well-defined, correct minimum of the likelihood function re- 
quires slow deweighting. 

2) Small tracking errors when aircraft parameters change re- 
quires fast deweighting. 

The first requirement exists because likelihood functions do not neces- 
sarily have minima in the correct place (i. e. , at c = c.) for short data 

2 - -t 

samples. It should be noted that this phenomenon is not inconsistent with 

maximum likelihood theory. The theory guarantees correct answers only 

asymptomatically. 
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The second requirement is less esoteric. As the aircraft changes flight 
condition, its parameters change in ramp-like fashion. The accumulated 
likelihood functions will then be out of date and cause parameter estimates 
to lag behind the true parameter values. The faster we deweight accumulated 
data, the less the lag. 

The design compromise for these requirements resulted in a time con- 
stant of 5 seconds. 

In addition to low-frequency deweighting of accumulated likelihoods, it 
was also found desirable to prevent very-high-frequency data (such as sharp- 
edged gust or command responses) from being accumulated. This was done 
by adding a relatively high-frequency, lowpass filter to the accumulation and 
highpass network. The total likelihood filter then takes the (digital) form 
shown in Figure 19, The symbols Ae^, Ae^, and denote generic inputs, 
intermediate states, and outputs, respectively. The filter is actually used 
to generate all likelihood functions, i = 1, 2, . . . N, and all components of 
VL and V 2 L. 


Ae. 



Unity gain lowpass 


Accumulation and 
highpass 


r = 0.6 s 
M 2 = ex P (-A/0.6) 

Figure 19. -Filter for likelihood functions. 


r = 5 s 
M = exp (-A/5) 
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Gain adjustment. -Adaptive control techniques that use explicit parameter 
identification require a gain schedule based on the parameter estimates. This 
section describes the gain adjustment portion of the PCMLE adaptive algorithm. 

Derivation of gains: Since the VALT performance specification was given 
in terms of desired command responses, it was straightforward to specify 
gains to match a desired closed-loop model. A gain matrix was computed at 
each flight condition based on the desired closed-loop characteristics. 

Plots of these gains versus flight condition (forward velocity) are contained 
in Appendix A (pitch axis) and Appendix B (lateral axis). Examination of the 
pitch-axis plots shows a strong dependence on forward airspeed for the 
majority of the elements. Simple functions of v (normalized velocity) were 
used to approximate the desired set of gains and thus define a "schedule. " 

The schedule was given in Table 4 and was checked at a number of flight 
conditions. 

The benefit of extending the gain calculation to a function of four variables 
was examined using a least-squares measure of fit. A "fit error" for each 
gain element was computed as the normalized sum of squares between the 
schedule and the desired gain value. The sum was done Over 33 flight con- 
ditions and the fit error was normalized by the fit error for zero gains. 

Results are given in Table 8. The first column shows the fit of the v 
schedule. (Note J = 1 for and K^, since the schedule is zero for V x 
feedback. ) Column 2 shows the improvement if all four of the parameters 
estimated by PCMLE are used for updating gains. Gain elements ^T5’ 

and Kgg shows some relative improvement. Columns 3, 4, and 5 of the table 
show the improvement in the fit if each of the three extra parameters are 
used one at a time with v. Consideration of these data shows that the improve- 
ment in fit is due to including parameter c^ (Sb^); thus, three gains will fit 
better if the basic v schedule is augmented with Cg dependence. The new 
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schedule was computed as: 


*12 = K 12 (v>- 0.6071 c 3 
*15 - *15 (v) + 0.582 Cg 
*26 * *26 (v) + 2-574c 3 


Next, the command responses with and without the Cg function were evaluated 
at a number of flight conditions. Evaluation of these time histories did not 
show any benefit of using Cg, Thus, although Cg improves the fit of several 
gains, this effect is not important in the closed- loop responses. Therefore, 
only the v dependence is retained. 


TABLE 8. -FIT ERRORS IN PITCH-AXIS GAIN SCHEDULE 


Gain 

J(v) 

J(v + 8) 

J(v + 6i) 

J(v + 62) 

J(v + 63) 

K(l, 1) 

1. o 

0. 61 

0. 88 

0. 75 

1. 0 

(1,2) 

0. 25 

0. 017 

0. 25 

0. 25 

0. 062 

(1,3) 

0. 033 

0. 0096 

0. 021 

0. 022 

0.029 

(1,4) 

0. 022 

0. 004 

0. 022 

0. 004 

0.020 

(1,5) 

0. 45 

0. 037 

0.45 

0.45 

0. 10 

(1,6) 

0. 020 

0. 002 

0. 02 

0. 002 

0. 02 

(2, 1) 

1. 0 

0. 93 

0. 98 

0. 98 

0. 98 

(2,2) 

0. 012 

0. 004 

0. 01 

0. 005 

0. 012 

(2,3) 

0. 83 

0. 70 

0. 72 

0. 83 

0. 83 

(2,4) 

0. 0155 

0. 0033 

0. 0158 

0. 0167 

0. 004 

(2,5) 

0. Oil 

0. 008 

0. 01 

0. 009 

0. 01 

(2,6) 

0. 37 

0. 20 

0. 37 

0. 36 

0. 27 


Note: PCMLE estimates v and 6, where 


6 


fia 33 

6 b 3i 

L 6b 32 
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Gain schedules as a function of parameter accuracy: In the VALT de- 
sign, gain adjustment is used to improve command responses and minimize 
cross coupling, not to maintain stability. Thus, it is logical to evaluate a 
"best" fixed-gain system. The best fixed gains can be determined by formally 
solving an optimization problem. 


The function being minimized by a weighted least- squares fit was ex- 
pressed as: 


= XL w.. (A + BK - A ).. 
jk jk m'jk 


where 

A, B = Vehicle model matrices 

A = Desired closed- loop matrix 
m r 

w^k = Weight factor for jk element of matrix (A - BK - AJ 
K = Gain matrix 


The matrix K can be determined by following the usual least- squares pro- 

p T 

cedure: set = 0 and solve for K. If only one flight condition is selected, 
this algorithm returns the set of gains for this condition. Summing over 
multiple flight conditions produces a weighted least-squares fit. In the fol- 
lowing analysis, the V =0 flight conditions were given three times the 

weight over V ^ 0 conditions in determining the gains, 
z 

The weighting factors were computed as follows. An individual weight- 
ing was applied to each of 16 gain elements. The weighting of each gain ele- 
ment was determined from the sensitivity of four responses to various errors 
in the closed- loop matrix. The four responses are: 
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r 2 

r 3 

r 4 


V due to V command (direct) 

z z 

0 due to V command (cross coupling) 
z 

V due to 0 command (cross coupling) 
z 

0 due to 0 command (direct) 


Define: 


6 


ij k 


Sri 



= Sensitivity of response i to change in 
gain element jk 

Evaluate at 
peak of re- 
sponse 


and 




(sum over the four responses) 


Then W, the matrix of elements w. k , is the weighting used in our least 
squares algorithm. 


It turned out that each a^ was important to only one of the four responses. 
Recall also that feedback is affecting only the V (row 2) and the q (row 3) 
rows of the closed-loop system matrix. The elements of the system matrix 
and their respective weighting are given below: 


Closed- loop parameters 


a 22 a 23 a 24 a 25 a 26 b 21 b 22 
a 32 a 33 a 34 a 35 a 36 b 31 b 32 


W 


0.04 0.25E-04 0. 25E-04 0.01 0. 6E-05 
16.0 0.01 0.01 4.0 0.25E-02 


0.06 0.25E-04 
144. 0.06 


54 



The elements of the W matrix are themselves interesting. The element, 
which is the feedforward from the 0 command to the collective and differen- 
tial collective, is the most sensitive. Errors in this parameter show up as a 
large AV z response to a 0 command. 

The gain schedules derived earlier in this section implicitly assume 
perfect knowledge of the scheduling parameter, v. If the uncertainty level 
of this parameter is known, it may be desirable to modify the gain schedule. 
(In an extreme case, the fixed gain configuration may perform better than 
a schedule if the scheduling parameter is in error.) This latter possibility 
was examined by computing the "best" fixed gains over an interval. The 
interval represents a ±2-sigma bound on the estimate. In the VALT design 
it was found that the gain schedule should be a function of the uncertainty level 
of the scheduling parameter. 

The accuracy of the scheduling parameter has been included by deter- 
mining a family of gain schedules for different accuracies of the parameter 
estimates. 

Three sets of gains were computed with the above procedure: 

• Low -sigma - Fit gains for flight conditions over an interval 

±0. 25 v about the nominal v. 

• High- sigma - Fit gains for flight conditions over an interval 

±0. 5 v about the nominal v. 

• Constant-gain - Fit gains over all flight conditions. 

Plots of the original v schedule, the functions for high- and low-v 
accuracy and the best constant value are contained in Appendix E. Note as 
the accuracy degrades (high sigma), the curves tend toward the best fixed- 
gain schedule. 
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Implementation: The preceding gain computation method has been imple- 
mented as a table lookup. The PCMLE v estimate and its accuracy estimate 
will be quantized and a gain matrix selected from the table. 

Performance on adaptive controllers on VALT simulation. -The PCMLE 
adaptive controller for the pitch axis was evaluated on the NASA-LRC VALT 
simulation. The performance agreed with linear simulations run at Honeywell 
in Minneapolis. This section contains a number of cases to illustrate the 
salient performance features of the algorithm. The PCMLE estimator uses 
a baseline set of three filters (channels) and estimates four parameters. 

The three channels Correspond to models at the following flight conditions: 

Channel 1 : v = 1. 0 

Channel 2: v = 0. 5 

Channel 3: v = 0 (hover) 

The PCMLE cases are summarized in Table 9. In these cases, the 
estimator was evaluated open- loop. A measured v was used by the gain 
schedule. 


TABLE 9. -VALT SIMULATION TEST CASES 


Flight condition 

Input 

PCMLE 
channel used 

Time history- 
trace 

Hover 

AB, step commands 

1.2 

Figs. 20, 21 

V = 19,5 m/s 

AB step commands 

3 

Figs. 22, 23 

X (63 ft/s) 



V = 41. 2 m/s 

AB , AV step commands 

1,3 

Figs. 24, 25 

(135 ft/s) 

z 


V v = 41. 2 m /s 

AB, gusts, sensor noise 

2 

Figs. 26, 27 

(135 ft/s) 



V = 79. 3 m/s 
(260 ft /s) 

AB, AV step commands 

2 

Figs. 28, 29 
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Simulation traces for each of these cases are shown in Figures 20 through 
29. Figures 20, 22, 24, 26, and 28 show the helicopter variables: 

0 - pitch attitude 

6g - differential collective command 

d' - collective command 
c 

- differential collective position 
q - pitch rate 

5 C - collective position 
x - forward velocity (earth reference) 

z - vertical velocity (earth reference) 

Figures 21, 23, 25, 27, and 29 show selected PCMLE variables: 

CH 

A 

V z 

A 

q 

T J(l) 

T J(2) 

T J(3) 

A 

Z 1 
A 

z 2 

Figure 20 shows that 6 and V command responses at hover meet speci- 

z 

fications. The step responses satisfy our criteria and there is no noticeable 
cross-axis response. Figure 21 shows that the likelihood function is mini- 
mum for channel 3 (hover). The * estimate should be zero and it is close 
to zero when estimating from either channel 1 or 2. This confirms the 
gradient calculations. 


channel (1, 2, or 3) being used for estimate 

estimate of V from selected channel 
z 

estimate of q from selected channel 

- | Likelihood functions for the three channels 

estimate of v parameter 
estimate of Sbgg parameter 
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Figure 22.-PCMLE time history test case 2 
helicopter variables. 
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Figure 23.-PCMLE time history test case 2- 
PCMLE variables. 
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Figure 24. -PCMLE time history test case 3 
helicopter variables. 
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Figure 25.-PCMLE time history test case 3 — 
PCMLE variables. 
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.31 m/s (1 ft/s) vertical gusts plus sensor noise 
PCM LE 

Trim = 41.2 (135 ft/s) 


Figure 26. -PCMLE time history test case 4- 
helicopter variables. 


64 



Figure 27. -PCMLE time history test case 4- 
PCMLE variables. 
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Figure 28. -PCMLE time history test case 5 
helicopter variables. 
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Figure 29. -PCMLE time history test case 5 
PCMLE variables. 
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Figures 22 and 23 show 0 responses at v = 0. 23 (37 kt). Channel 3 
still fits best as expected. The v estimate is a little low, but still acceptable. 

Figure 24 shows 6 and V command at a v = 0. 5 (80 let). Again the com- 

z 

mand responses are good and there is no cross-axis response. Figure 25 
shows v estimated from channel 1 and channel 3. The parameter v should 
be 0. 5 at this condition and the estimate is reasonably close. Channel 2 
should have the minimum TJ, but as a result of errors in other parameters, 
TJ(1) and TJ(2) are very close. 

Figures 26 and 27 present a repeat of the v = 0. 5 flight condition with 
gusts and sensor noise. The v estimate is slightly lower than 0. 5, but accept- 
able. Channel 2 should be minimum, but the likelihood function of channel 
3 is slightly smaller. 

Finally, Figures 28 and 29 show the performance at the maximum 

velocity, v = 1 (160 kt). Again, 6 and V responses are very good, with a 

z 

barely detectable cross -axis response. The estimate of v is close to 1 as 
it should be. Again, channel 3 is minimum, although channel 1 should be at 
this condition. 

This section has presented a design based on explicitly identifying four 
CH-47 pitch-axis parameters (v plus perturbation parameters on the B(3, 1), 
B(3, 2) and A (3, 3) coefficients). This section described the analysis used to 
select the parameters to be identified and the channel structure to be used 
in the algorithm. Next, the use of the v estimate in scheduling the gains was 
considered. The resulting gain schedule was implemented as a table lookup 
that uses the v estimate as well as its indicated accuracy. Finally, the per- 
formance of the identifier on the VALT simulation was considered. The 
accuracy of identifying v was consistent with predictions. It was observed 
that in some situations the minimum channel is not at the correct value of 
v (even though the estimate is reasonable). This result is investigated in 
Section 10 where the results of processing actual flight data are presented. 
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SECTION 7 

MODEL REFERENCE ADAPTIVE ALGORITHM 


The block diagram of the model reference (MR) algorithm as applied 
to the CH-47 longitudinal axis is shown in Figure 30. Note the command 
model is lower order than the vehicle model. 


Overview of Algorithm 


Consider the CH-47 longitudinal axis in the form given in Section 4: 
x = Ax + B6 


and a model 


x = A x + B u 
m m m m 


where 


X = <V x , V 2 , q, 0) 

T 

5 = (6-n > 6 ) are the actuator commands 

Jl> c 

T 

u = (V , 0 ) are the command inputs from guidance 
z c 
c 


X = (V , q , 0 ) 

m z m m 

m 


The three-state model provides the desired 0 and V command 

o z 

c 

responses. The 0 response is second-order, the V response is first- 

z m 

order. Including integral control on each command error adds two states 
to the control law. 
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V z q 9 1 

Note: j K^], K x is partitioned to s« 

i 

Figure 30. -Structure of MR a 
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The control. Law has the form discussed in Section 5: 


6 = K x x + K u u 

which gives a closed- loop: 

x = (A + BK ) x + BK u. 

X LI 

Here x includes the two command error integrator states. 

The function of the adaptive algorithm is to adjust the gains so the 
closed-loop matrices match the model. 

As A and B vary with flight conditions, K and K need to be continually 

X LI 

adjusted for the plant to match the model. The gain matrices K x and K 
contain a total of 10 individual terms. 


In the MR algorithm, the gains are adjusted to force a Liapanov func- 
tion of model-following errors and parameter errors to zero. It is well- 

known that a Liapanov-derived gain adjustment rule is stable in the error, 

3 

e. If the parameter errors do not go to zero, they remain bounded. 

The first form of a Liapanov function used in the study was: 

V = 1/2 e T Pe + 1/2 Tr(AAAA T ) + 1/2 TR(ABAB T ) 
where A A and AB are the difference between the model and the vehicle. 


Taking the appropriate time derivatives yields adjustment rules of the 


form; 

* T 

A A = -Pex 1 

. rp 

AB = -Peu 1 
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where x is the measured state, u is the input and e is the error between 
model and plant states. 


Since the adptation is implemented as an adjustment of feedforward 
and feedback gains, the adjustment rules need to be expressed as AK func- 
tions rather than direct changes to the plants A and B matrices: 


^(A+BK x ) = -Pe T x (3) 

K = -B" 1 Pe T x 

a? (BK u> ■ - peTu 

* 1 T 

K u = -B Pe u (4) 


The above form implies a B"^ can be found that works sufficiently 
well over the flight envelope, since in this application B also varies with 
flight conditions. 


The P matrix is determined from usual stability considerations: 


m 


P + PA = Q * 0 


m 


(5) 


If Q = \Q q , where \ is an arbitrary scaling parameter, then the P 
matrix in Equations (3) and (4) is also scaled. Thus, P is a somewhat 
arbitrary design parameter [as long as Equation (5) is satisfied] . 


In past applications, problems with the gain adjustment equations arose 

2 

because they are sensitive to the signal levels of e, x, and u. (If a P was 
found that worked well with large commands, the adaptation was too slow 
with small inputs; if tuned to work well with small inputs, then it became 
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very oscillatory with large inputs and usually diverged due to numerical 
difficulties.) 

It was recognized that monotonic functions of the basic Liapunov func- 
tion could also be used and these forms could be used to "normalize" signal 
levels in the K equations. Thus, by using different Liapunov functions, 
different forms of the gain adjustment equations are obtained. The three 
forms investigated are given in Table 10. The third form is the one finally 
selected for implementation. 


TABLE 10. -MODIFIED LIAPUNOV FUNCTIONS 


Liapunov function . 

Adjustment rule 

V = ePe T + 1/2 Tr AA AA T 
+ 1/2 TrABAB 

• — 1 rri 

K x = -B Pe x 
K 2 = -B -1 Pe T u 

V = 1/2 In (e T Pe + 1) 

+ 1/2 TrAAAA T + 1/2 TrABAB T 

K nT Pe T x 

1 (e i Pe + 1) 

k„ Pe T u 

£ (e i Pe + l) 

V = ^ T Pe + 1 

+ 1/2 TrAAAA T + 1/2 TrABAB T 

K. = - — — * Pe T x 

\^ T Pe + 1 

k. - - B' 1 Pe T u 

\^ T Pe + 1 
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Simulation Results 


A number of experiments were conducted using the CH-47 longitudinal 
model to develop the MR algorithm. This included considerations for 
model-following errors caused by states not being controlled and the ef- 
fects of various design parameter on performance. 

The model is 



2 

s + 2 



7 

s 2 + 4s + 7 


w = 7 rad/s 
p = 0.76 


P = \ 


0.25 0 0 

0 0. 1944 0. 0833 

0 0.0833 1. 4167 


\ = 300 


A fixed set of and Kg is used as starting gains at each flight con- 
dition. (They are given in Figure 30. ) These gains then adjust based on 
the update rules. and Kg are presently adjusted with the rule shown in 
row 3 of Table 11. 
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TABLE 11. -VALT SIMULATION TEST CASES 


Flight 

condition 

Input 

Model 

reference 

mode 

Time 

history 

trace 





Hover 

V z , B step commands 

Fixed gain 

Fig. 31 

V = 41. 2 m/s 
x (135 ft/s) 

V , Q step commands 

Fixed gain 

Figs. 32, 33 

V = 79. 3 m/s 
x (260 ft/s) 

V , 0 step commands 

Fixed gain 

Figs. 34, 35 

V = 79. 3 m/s 
x (260 ft/s) 

G $ V z step commands 

Adapting 

Figs. 36, 37 

V v = 79. 3 m/s 
(260 ft/s) 

0, Gusts and sensor 
noise 

Adapting 

Figs. 38, 39 


Performance of MR Algorithm on VALT Simulation 

The performance of this algorithm on NASA-LRC's simulation is sum- 
marized in this subsection. This algorithm adjusts gains based on errors 

between the model states and three measured helicopter states (V , q, 0). 

z 

The starting set of gains was selected to work well at hover. In practice, 
they work well up to forward speeds of 41. 2 m/s (80 kt). The major perfor- 
mance degradation occurs at the maximum forward speed of 79. 3 m/s (160 
kt) and is manifested primarily by large cross-axis responses. Thus, the 
model reference adaptation is shown at the high-speed condition, since its 
effects are most apparent there. 

Simulation traces for these cases are shown in Figures 31 through 37. 
Figures 31, 32, 34, and 36 show the same helicopter variables as the PCMLE 
runs do. Figures 33, 35, and 37 show the following variables: 
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35 (20) 


0 

rad 0 

(deg) 


.79 (45) 


0 

rad 0 

(deg) 


3.41 (180) 


* 

rad 0 

(deg) 


4 


In. 


0 


.17 (10) 


q 

rad/s 0 

(deg/s) 


5 C 

in. 


4 

0 


v 31 (100) 
x 

m/s 0 

(ft/s) 


6.2 ( 20 ) 


V z 

m/s 0 

(ft/s) 



Trim =3.1 m/s (10 ft/s) 


Figure 31.- MR time history test case 1 
helicopter variables. 
















3.1 ( 10 ) 



Figure 37. -MR time history test case 4 — command input 
model state, and error variables. 


U - V command input to model 

v z 

LL. - 0 command input to model 

U 

V - V modeL state 

z z 

m 

0 - 0 model state 

m 

E - V error between model state and helicopter 
v z 

z 

E q - q error between model state and helicopter 

E - 0 error between model state and helicopter 

0 

It should be noted that the inputs to the model are high-passed so steady- 
states do not appear on the traces. (Basically, this is required to remove 
trim effects. ) 

Figure 31 shows V and 0 command responses at hover. They meet 

z 

requirements and show no cross-axis response. 

Figures 32 and 33 show 0 and V command responses at the 41. 2 m/s 

z 

(80 kt) and 79. 3 m/s (160 kt) with the starting gains fixed. At the 79. 3 m/s 
condition, some cross-axis response is present; otherwise, the command 
responses are reasonable. 

Figures 34 and 35 show the 79.3 m/s condition with the gains adapting. 
Note the reduction in the cross-axis response for 0 command. Also note 
that the V z command that follows has a very slow response because some 
gain values moved to accommodate the 0 command. 

Figures 36 and 37 illustrate the fundamental problem with this algo- 
rithm. The system is at the high-speed (79. 3 m/s) condition with gusts 
and sensor noise. After operating for several minutes, a sequence of 0 
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commands is applied. The random disturbances had caused the gains to 
drift sufficiently that an unstable combination resulted. 

This section has presented the design of a model reference algorithm. 
The structure of the control law is identical to the adaptive design of Sec- 
tion 6. The difference lies in the method used for determining the gain 
matrices. In this design, the gain is adjusted with a Liapunov procedure 
to drive the response errors toward zero. Simulation studies indicated 
good performance except when measurement noise is present. In this situa- 
tion, the gains continue to drift in response to the measurement noise. 

This result is discussed further in the next section which compares the two 
candidate adaptive designs. 
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SECTION 8 

CONCEPT COMPARISON AND SELECTION 


The two designs from Sections 6 and 7 are compared in Table 12. This 
comparison summarizes the simulation performance, the design process, 
the complexity of the algorithm, and its simulation performance. 

The PCMLE adaptive design was recommended for flight evaluation 
based on its superior performance with gust and sensor noise present. In 
addition, the digital formulation of the algorithm offers better numerical 
stability than discretization of the model-following equations. 

The PCMLE algorithm for real-time parameter estimation is relatively 
complex compared with current flight control software. The various PCMLE 
subroutines are written in FORTRAN with an emphasis on speed of execution. 
During development, the software was run on a CDC-6600 and a Honeywell 
H-6080 computer in Minneapolis. 

Load maps from CDC-6600 compilations show that about 5K of memory 
are required (about 1. 6K of this represents labeled common blocks). Timing 
checks made on the H-6080 indicate a maximum of 5.6 ms, with minor 
variation between the five subcycles. 
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TABLE 12. -concept comparison 


Parameter 

PCMLE adaptive 
controller 

Model reference 
adaptive controller 

a) Command responses 

(av z , A e) 

Both designs satisfy response specifications. Command 
responses are nearly identical for nominal conditions 
(no gusts or sensor noise). 

b) Stability of control- 
ler 

Equal for nominal operation. 

c) Stability of adaptive 
portion of algorithm 

Estimator stable and 
unbiased. Only plants 
that do not exhibit 
approximately linear per- 
turbation responses would 
cause problems. 

Theoretical. In practice, 
numerical problems can 
cause algorithm to diverge. 
Related to item f below. 

d) Effects of turbulence 

No problem. It is accounted 
for in Kalman filter design. 

Causes drift in some gain 
elements. 

e) Effects of sensor 
noise 

No problem. It is accounted 
for in Kalman filter design. 

Causes drift and biases 
in some gain elements. 

f) Implementation 
philosophy 

a) Explicit gain adjustment. 
Requires designing a 
gain schedule based on 
parameter estimates. 

b) Direct digital formulation. 

c) Offers the lowest sample 
rate. 

a) Implicit gain adjustment. 

b) Discretization of analog 
equations. Problems in 
integrating continuous 
equations with large At. 

c) Problem not we 11- posed 
if perfect model-follow- 
ing is not possible. 

g) Complexity (real- 
time code only) 

1602 core FORTRAN code. 
1608 core common blocks. 

452 core FORTRAN code. 
256 core common blocks. 

h) Computation time 

6 ms per cycle (3 channels - 
4 parameters estimated). 

Less than 1 ms per cycle. 
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SECTION 9 

GUIDANCE INTERFACE 


As discussed previously, the controller was designed for 0, 0 , t// , and 

V commands. Since the existing VALT guidance algorithm generates 
z 

velocity commands, an interface was required to convert velocity commands 
to attitude commands. The guidance interface is shown in Figure 38. 


Lateral velocity 
error 




0 N 


s/0.5 + 1 
s/5.0 + 1 

K* = 0.02 

K, =0.2 
K e = -0.04 



Figure 3 8. -Guidance interface. 


The bank-angle command is generated from a sum of lateral velocity 
error and the nominal bank angle computed by the guidance law. 


The pitch- attitude command is more difficult to generate, since a 
nominal 0 is not available. The longitudinal velocity error is applied to a 
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proportional, plus integral and used as the attitude command. A lead- lag 
filter was added to improve the dynamic response of the coupled guidance/ 
control system. A root-locus analysis predicted the closed-loop response 
to a velocity command should have a damping ratio of 0. 5. 

A time history of trajectory-following using the attitude controller and 
the guidance coupler is shown in Figure 39. The performance is satisfactory. 
An overshoot of approximately 20 percent can be observed in the forward- 
velocity trace as the vehicle completes the trajectory and comes to a hover. 
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SECTION 10 

ESTIMATION OF PARAMETERS FROM FLIGHT DATA 


The PCMLE algorithm was evaluated against prerecorded flight data 
to assess the performance of the explicit identifier. Multiple iterations of 
the Newton-Raphson parameter steps were used to obtain the best possible 
convergence. This section describes the sensor data, the procedure used 
to obtain parameter estimates, and draws conclusions from.the results. 

Analysis of Sensor Data 

Recorded flight data for a number of maneuvers was available from an 
earlier flight test program. The variables of interest for the PCMLE algo- 
rithm are vertical velocity, pitch rate, pitch attitude, and the two controls, 
differential collective and collective. A functional diagram of the pitch 
axis is shown in Figure 40. The variables Z13, Z16, and Z19 were used 
to obtain approximations for 6 C and 6^ by assuming the actuators have unity 
transfer functions. 

Thus, 6 = Z19 and 6-c = Z16 + Z13. 

C -D 

The flight data, which had been recorded at 40 sps, was highpass fil- 
tered to eliminate unwanted trim values in both control inputs and measure- 
ments. (This highpass operation is part of the PCMLE algorithm.) Initial 
processing was performed to validate scale factors and establish reasonable 
sensor noise levels. 

Since the stability augmentation system (SAS) had not been cancelled in 
these flight tests, its effect had to be included. The recorded variable Z13 
was compared with an estimate computed from the SAS transfer function and 
the measured pitch rate. 
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V z 

0 

q 


Figure 40. - Recorded sensor outputs from flight records. 


The two signals differed by a scale factor which changed with maneuver. 
The decision to retain the recorded version was made following test cases 
in which its use resulted in better matches to the recorded measurements. 
Further analysis of the measurement residuals showed that in these calibra- 
tion runs they were dominated by the vertical velocity measurement. Re- 
moval of this term from the observation equations resulted in acceptable 
residuals for the two remaining sensors: pitch rate and pitch attitude. The 

estimation process was insensitive to assumed noise statistics for these 
measurements over a wide range of values. 

A processing rate of 10 sps (every fourth data point) was investigated 
to reduce computation time. Although this rate seemed well within required 
sampling limits, there were significant residuals that indicated that a higher 
rate was required. Therefore, the original sample rate of 40 sps was used. 
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Procedure 


The basic MLE algorithm used in the identifiability analysis was used. 
Some numerical problems were encountered in finding solutions to the 
Ricatti equation. This occurred since our algorithm for solving the Ricatti 
equation requires that the initial set of Kalman gains stabilize the system. 
This was often difficult to do in view of the parameter adjustments occurring 
at each iteration. In retrospect, it would have been more efficient to solve 
the Ricatti equations with a method such as Potter's algorithm which removes 
the requirement for initial stabilizing gains. 

Because of the apparent small gust levels, a number of runs were made 
using zero filter gains. Using zero gains often caused the parameter errors 
to diverge. This occurs because the pitch-axis model is unstable at a num- 
ber of flight conditions. In this situation, the residuals diverge and the 
likelihood functions have such steep sides that it is very difficult from a 
numerical viewpoint to locate the minimum with Newton-Raphson steps. 

The "zero gain" procedure was modified as follows. 

Once an initial set of parameters was obtained, gains were added, with 
the hope that parameter updates would be small enough to allow use of the 
same stabilizing initial gains. In most cases this proved to be true. How- 
ever, for the final parameter identification experiments, the "zero-gain" 
procedure was not used. 

The ill-behaved nature of the likelihood function necessitated careful 
attention to the allowable search directions in maximization. After a number 
of unsuccessful attempts at parameter localization with various data runs, 
it was decided to abandon v" (forward speed estimate) as the fundamental 
parameter for purposes of function optimization. Perturbation parameters 
on a^g and b^ appeared to be much more closely related to the quality of 
fit. In the final set of parameter identification runs, v was set to nominal 
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value and a 33 , t> 31 were allowed to find their best values. Once localized 
in the parameter space in this manner, all the other perturbation parameters 
plus initial conditions were estimated. The final model used in the Kalman 
filter consisted of a fourth-order state representation of the longitudinal- 
axis dynamics with two control inputs as in the baseline PCMLE algorithm. 

A third "control" was added to simulate the effect of highpass filtering the 
contribution due to the initial conditions. 

The final noise model was developed to account for uncertainties in the 
control inputs and the two sensors used (pitch rate and pitch attitude). State 
noise with an rms value of 2 percent of the control matrix was added to 
account for random perturbations of the recorded control inputs. Response 
measurements of pitch rate and pitch attitude were assigned noise statistics 
CTq = 0.01 rad/s and 06 = 0.01 rad, respectively. 

Results 

The final set of 12 maneuvers analyzed is defined in Table 13. 

The procedure described resulted in the parameter estimates shown 
in Table 14 for each of the 12 maneuvers analyzed. The quality of fit is 
rather good. Plots of fit for each maneuver are contained in Appendix F. 
Three plots are used to describe each maneuver. The first plot shows the 
two control inputs. The second and third plots show the highpassed measure- 
ments of 0 and q and their corresponding residuals for the final iteration. 

It should be emphasized that only two parameters significantly improve over 
their a priori uncertainty. Table 15 illustrates the reduction in parameter 
uncertainty for a typical case (maneuver 7). Note that only a 33 and b 31 
(plus two initial conditions) are significantly improved. 
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TABLE 13. -MANEUVERS ANALYZED WITH 
FLIGHT DATA 


No. 

Data run 

Serial no. 
of file 

Flight condition 

V x (m/s) 

V z (m/s) 

1 

5 

34007 

0 

0 

2 

9 

31017 

0 

0 

3 

13 

34031 

0 

0 

4 

17 

35001 

10. 3 

+7. 63 

5 

25 

31005 

20. 6 

0 

6 

27 

35008 

20. 6 

-5. 08 

7 

29 

31017 

20. 6 

-2. 54 

8 

31 

31023 

20. 6 

+2. 54 

9 

33 

31026 

20. 6 

+5. 08 

10 

39 

33010 

30. 9 

0 

11 

41 

30010 

30. 9 

+2. 54 

12 

43 

33012 

30. 9 

+2.54 


The basic conclusions are that near the correct value of v, good param- 
eter estimates can be obtained. If v is wrong, the algorithm may find 
another local minimum and erroneous parameter estimates are obtained. 

It appears that reasonably good fits can be obtained at a number of values of 
v. When the baseline PCMLE algorithm was simulated, similar results 
were obtained. 

Based on our analysis of flight data, we have concluded that v is pro- 
bably not a good dominant parameter. If the baseline PCMLE algorithm 
had been flight-tested, it probably would have experienced problems with 
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TABLE 14. -parameter values for pitch -axis identification 


No. 

V 

a ll 

a 12 

a !3 

a 14 

a 2l 


a 23 

a 24 

a 31 

a 32 

a 33 

a 34 

b ll 

b 12 

b 21 

b 22 

b 3 1 

b 32 

1 

0 

-0. 0181 

-0.000134 

2.41 

-32.2 

-0. 180 

-0.499 

0. 00594 

-0. 0685 

0.0113 

0. 000324 

0.460 

0. 981 

0. 121 

-0. 000818 

0. 0000724 

-7.80 

0.272 

0. 000119 

2 

0 

-0. 0181 

-0.000069 

2. 41 

-32. 2 

-0. 180 

-0.500 

0. 00151 

-0. 0141 

0.0144 

0.0000988 

0.313 

0. 974 

0. 121 

0 

-0. 0000254 

-7.80 

0.272 

0 

3 

0 

-0. 0181 

-0.000136 

2.40 

-32.2 

-0. 180 

-0. 499 

0. 00281 

-0. 0293 

0. 0109 

0.000308 

0. 369 

0. 929 

0. 120 

-0. 000562 

0.0000347 

-7.80 

0. 277 

-0. 0115 

4 3 

0. 125 

-0. 0186 

-0.000814 

2.41 

-32.0 

-0. 150 

-0. 494 

0. 00231 

-22. 4 

0.000823 

0. 00227 

-0. 167 

1. 46 

0. 120 

-0.000703 

0.0870 

-8. 18 

0. 269 

0. 0289 

5 

0. 25 

-0. 0202 

-0.00318 

2. 40 

-31.3 

-0. 120 

-0. 549 

0. 00108 

-45. 3 

0. 00764 

-0. 00333 

0. 128 

1.07 

0. 121 

0. 00953 

0. 176 

-8.60 

0. 285 

-0. 0110 

6 

0,25 

-0. 0198 

-0. 00403 

2. 42 

-31. 3 

-0. 122 

-0. 825 

-0. 0144 

-53,7 

0. 0216 

0. 0000202 

0. 198 

2. 12 

0. 119 

-0. 163 

0.211 

-8. 20 

0.298 

0. 0189 

7 

0. 25 

-0. 0203 

-0. 00328 

2.42 

-31. 3 

-0. 120 

-0.498 

0. 00201 

-45. 0 

0. 0201 

0.000373 

0.284 

2. 18 

0. 122 

0. 00228 

0. 175 

-8.55 

0.295 

0. 00597 

8 

0. 25 

-0. 0203 

-0.00368 

2. 40 

-31. 3 

-0. 120 

-0. 989 

-0. 00705 

-43.7 

0.0142 

-0.0186 

0. 151 

1. 19 

0. 124 

-0.0189 

0. 136 

-0.40 

0. 265 

-0. 00277 

9 

0. 25 

-0.0203 

-0.00429 

2. 40 

-31. 3 

-0. 120 

-1. 19 

-0. 0256 

-42. 7 

0.0222 

-0.0245 

0.209 

0.783 

0. 128 

-0.0164 

0. 117 

-9.83 

0.258 

0. 00740 

io a 

0. 375 

-0. 0229 

-0.00286 

2. 37 

-30. 2 

-0.0899 

-0. 0671 

-0.0664 

-66. 4 

-0. 00499 

-0.00114 

0.282 

2. 56 

0. 121 

0. 158 

0.242 

-8. 84 

0. 295 

0. 0409 

ll a 

0. 375 

-0. 0229 

-0.00286 

2. 37 

-30. 2 

-0. 0899 

-0. 0671 

-0. 0664 

-66.4 

-0. 00499 

-0.00114 

0.282 

2.56 

0. 121 

0. 158 

0.242 

1 

CO 

CO 

►u 

0.295 

0. 0409 

12 a 

0. 375 

-0. 0229 

-0.00286 

2. 37 

-30.2 

-0. 0899 

-0.0671 

-0. 0664 

-66. 4 

-0.00499 

-0.00114 

0.282 

2.56 

0. 121 

0. 158 

0.242 

-8.84 

0.295 

0. 0409 


a Unable to find filter gain solution for next update. Error has probably not been minimized. 


CD 



Parameter 

CT o 

Of 

a n 

0. 0045 

0. 00460 

a 12 

0. 009 

0. 00914 

a 13 

0. 4 

0.401 

a, . 

3.0 

3.00 

14 



a 21 

0.04 

0. 0400 

a 22 

0.25 

0. 250 

a 23 

1. 6 

1. 600 

a 24 

10.0 

10. 0 

a 31 

0. 02 

0. 0531 

a 32 

0. 006 

0.00913 

a 33 

1. 25 

0. 102 

a 34 

0. 5 

1. 82 

b ll 

0. 06 

0.0609 

b 12 

0. 3 

0. 322 

b 21 

0.25 

0. 250 

b 22 

0. 8 

0. 800 

b 3 1 

0. 25 

0. 0471 

b 32 

0. 08 

0. 119 

X0(1) 

1. 0 

0. 855 

X0(2) 

1. 0 

1. 00 

X0(3) 

0.017 

0. 133 

XO (4) 

0.017 

0. 00873 







local, minimum as our off-line processing did. Perhaps, better parameter- 
ization would use the control effectiveness parameters (bgp b,^) as dom- 
inant, since they are most easily identified. 

One of the reasons v was selected as the dominant parameter was that 
analysis of the gains over the flight envelope showed a v dependence. In 
fact, a gain schedule based on velocity was evaluated and it performed well. 
The performance differences of the attitude control system with and without 
gain scheduling are primarily in the area of decoupling. Use of a gain sched- 
ule improves the decoupling of the A0 and AV responses. Since the de- 
coupling is dominated by the control input matrix, it is expected that a gain 
schedule could be expressed in terms of the B matrix. This would then be 
a more natural schedule to use with an explicit parameter identifier. 


97 


SECTION 11 
CONCLUSIONS 


Good command responses exhibiting desired speed of response and 
minimal cross coupling were obtained in the pitch axis with simple gain 
schedules based on forward airspeed. In the lateral- directions l axes no 
gain adjustment was required. These control laws were successfully flight 
tested at the Wallops Island facility using the CH-47. 

Use of velocity guidance with an attitude command system generally 
did not result in as tight a trajectory-following capability as maybe desired. 
The guidance algorithm should be designed in conjunction with the control 
system for optimum performance. 

Adaptation for vehicles like the CH-47 is marginally useful and mar- 
ginally successful. The identifiability analyses based on computed Cramer- 
Rao lower bounds showed that only about 7 of 26 longitudinal parameters 
can be extracted from the sensor data. Relatively large maneuvers are 
required to obtain significant improvement over a priori knowledge. The 
dominant scheduling parameter (forward airspeed) can be found to an ac- 
curacy no better than 20 to 30 kt (about 15 to 20 percent of maximum speed). 
Off-line evaluation of the identifier supported these identifiability results. 
This evaluation also indicated that the fundamental parameter identification 
problem is plagued with local minima. Initial estimates must be reasonably 
close in order to converge to the correct parameter set. 

These experiments suggest real-time parameter estimation is feasible 
in modern flight control. Whether it offers benefits for a particular applica- 
tion depends on the nature and magnitude of plant uncertainties and alternate 
methods available to achieve successful control. 
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APPENDIX A 

GAINS FOR LONGITUDINAL AXIS 


Actuator commands have the form: 

6 = K x + K u 
xu 

This appendix contains design values for K and K for 33 flight con- 

X u w 

ditions plus the v gain schedule. The design values were computed using the 
procedure described in Section 5. The two columns of matrix K u equal col- 
umns 5 and 6 of matrix K so only K is plotted. Thus, individual terms are 
denoted as 



The data, in the form of computer plots, are presented in Figures 41 
through 52. 
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v, Normalized forward speed 


Figure 41. -Element Kjj versus v -- longitudinal axis. 
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Figure 45. -Element K 1& versus v -- longitudinal axis. 
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Figure 47. -Element K 1? versus v -- longitudinal axis. 
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Figure 52. -Element Kgg versus v 


longitudinal axis. 
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APPENDIX B 

GAINS FOR LATERAL-DIRECTIONAL AXES 


Actuator commands have the form: 


6 = K x x + K u u 

This appendix contains design values for K and K for 33 flight con- 
ditions plus the v gain schedule. The design values were computed using 
the procedure described in Section 5. The two columns of matrix K u equal 
columns 6 and 7 of matrix K so only K is plotted., Thus, individual terms 
are denoted as 



K 11 

CM 

H 

K 13 

K 14 

K 15 

K 16 

K 17 

K 21 

K 21 

K 23 

K 24 

K 25 

K 26 

K 27 




The data, in the form of computer plots, are presented in Figures 53 
through 66, 
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Figure 53. -Element K. , versus v 


lateral-directional axes 
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Figure 57. -Element versus v -- lateral-directional axes. 
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Figure 5 9. -Element K 1? versus v -- lateral-directional axes. 
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APPENDIX C 

LATERAL-DIRECTIONAL TIME HISTORIES 


This appendix contains time histories of the lateral-axis controller in 
response to 6 and t// commands. The traces were generated with the VALT 
simulation at NASA-LRC. Three flight conditions are given: hover, 41.2 
m/s forward velocity, and 79. 3 m/s forward velocity. 

The data are presented in Figures 67, 68, and 69. 































APPENDIX D 

CH-47 MODEL PARAMETERS 

This appendix contains plots of the elements of the A and B matrices 
for the longitudinal axis of the CH-47. The parameterization of these ele- 
ments in terms of v (normalized forward airspeed) is also shown. 

The vehicle model has the form (see ref. 4): 

a ll a 12 a 13 a 14 V x b ll b 12 

a 21 a 22 a 23 a 24 V z b 21 b 22 

+ 

a 31 a 32 a 33 a 34 q b 31 b 32 

0 0 10 e 

The data, in the form of computer plots, are presented in Figures 70 
through 87, 
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Figure 72. Element a 13 versus v -- CH-47 longitudinal axis. 
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PARAMETERIZATION 


Appendix D 
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Figure 75. -Element a 22 versus v -- CH-47 longitudinal axis. 
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Figure 76. -Element a 23 versus v -- CH-47 longitudinal axis. 
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Figure 78. -Element a g j versus v -- CH-47 longitudinal axis. 
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Figure 80. -Element a^g versus v -- CH-47 longitudinal axis. 
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Figure 82. 


-Element versus v 


CH-47 longitudinal axis. 
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Figure 83. -Element versus v 


CH-47 longitudinal axis. 


145 
















2 5000E-0 1 3.75000E-01 4.25000E-01 4.75000E-01 


Appendix D 











APPENDIX E 

LONGITUDINAL- AXIS GAIN SCHEDULE 
WITH ESTIMATED VELOCITY 


This appendix presents the ideal gain schedule, the best fixed gains, 
and the gain schedules for two levels of uncertainty in the scheduling param- 
eter (v). This analysis is described in Section 6. 

The data, in the form of computer plots, are presented in Figures 88 
through 101. 
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Figure 97. -Element K 24 versus v 


longitudinal-axis gain schedule. 
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Figure 100. -Element 


versus v 


longitudinal-axis gain schedule. 
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APPENDIX F 


PARAMETER ESTIMATION USING 
RECORDED FLIGHT DATA 


This appendix presents plots of 12 maneuvers analyzed in Section 10. 
For each maneuver the two controls, pitch rate and its residual and pitch 
attitude and its residual, are plotted. These plots have been normalized. 
The normalization factors are 

Controls (6_, 6 ): 1 = 2,54 cm 

D L 

Pitch rate: 1 = 1 rad/s 

Pitch attitude: 1 = 1 rad 

The data, in the form of computer plots, are presented in Figures 102 
through 137, 
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Figure 106, -Pitch rate versus time -- maneuver no. 2. 
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Time, s 

Figure 109. -Pitch rate versus time -- maneuver no. 3. 
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Control Input Pitch attitude 
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Figure 113. -Pitch attitude versus time -- maneuver no. 4. 



Time, s 


Figure 1 14. -Control input versus time -- maneuver no. 5. 


171 













Pitch attitude Pitch rate 


Appendix F 



Time, s 

Figure 121. -Pitch rate versus time -- maneuver no. 7. 



Time, s 

Figure 122. -Pitch attitude versus time -- maneuver no. 7. 
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Pitch attitude 
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Figure 129. -Control input versus time -- maneuver no. 10. 



Figure 130. -Pitch rate versus time -- maneuver no. 10. 
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Figure 133. -Pitch rate versus time -- maneuver no. 11. 
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Figure 134. -Pitch attitude versus time -- maneuver no. 11. 
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